Colección de ejercicios resueltos de 


Matemáticas Aplicadas a la Biología 


En estas páginas se expone una colección de ejercicios resueltos de la asignatura Matemáticas 
Aplicadas a la Biología. Todos ellos provienen de exámenes de los últimos cursos. 


De la misma manera que no se aprende a jugar al tenis mirando partidos por la televisión, 
no se aprende a resolver problemas leyendo problemas resueltos. Es necesario enfrentarse al 
problema, intentando resolverlo una y otra vez, si fuera preciso, acudiendo a la solución sólo 
para comprobar el resultado, o si no se consigue hacerlo. 


Las Matemáticas utilizan un lenguaje que necesariamente debe ser claro, conciso y, sobre todo, 
preciso. No consiste este lenguaje en un conjunto de fórmulas puestas unas al lado o debajo de 
las otras. Es necesario expresar, mediante símbolos o mediante frases, la lógica del razonamiento 
que se está realizando. 


Los resolución de los ejercicios de esta colección está escrita (en casi todos los casos) con todo 
detalle con un doble objetivo: 

(a) facilitar su comprensión; 

(b) mostrar cómo se debe escribir la solución de un ejercicio: de alguna manera es lo que 
los prefesores desearíamos encontrarnos al corregir un examen. 


Ejercicios del Tema 1: Repaso de 
resultados previos 


1. En un cultivo de laboratorio, el número de bacterias (medido en millones) durante las primeras 100 horas 
viene dado por 
N(t) =25+te 1%, te [0,100]. 


Determinar los periodos de crecimiento y decrecimiento de la población así como los momentos en que ésta 
alcanza su máximo y su mínimo absolutos. 


Comenzamos por calcular la derivada de la función N (t): 
-1 t 

N'D = e7t/ 44 sia ) 2 P E ) 
(t) =e + te To e — 


t 
que se anula si 1 — — = 0, es decir, para t = 10. 


71/10 es positivo Vt, es fácil ver que 


l N'(t) > 0 para t € (0,10) l N es creciente en (0, 10) 
=> 


Puesto que e 


N'(t) < 0 para t € (10, 100) N es decreciente en (10, 100) 


En consecuencia, N tiene un máximo local en t = 10, que también es máximo absoluto, ya que es monótona a 
ambos lados de este punto. 
El mínimo absoluto hay que buscarlo en los extremos del intervalo, al no haber mínimos locales: 


100 
N(0)=25 y  N(100)=25+ 100e7 100/10 — 25 + z ~ 25.0045 


Por lo tanto, el mínimo absoluto se produce en t = 0. 


2. El rendimiento de una determinada plantación de árboles viene dado por 


donde x es la distancia en metros entre los distintos árboles. ¿A qué distancia se deben plantar unos árboles 
de otros para conseguir una mayor producción? 


La distancia óptima será la que produzca el mayor rendimiento, es decir, la que produzca el máximo absoluto 
de la función f(x) en el intervalo (0, +00), ya que la distancia debe ser positiva. 


a?—8 
4 


Maximizar f(x) = 


para x € (0, +00) 
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f es continua y derivable en (0, +00). 
Derivada: 
Fa) = 2g- zt — 4r? (x? —8) 2g? — 4r? +32 32-22? 


g? g’ z5 
Puntos críticos de f: 


f'(£)=0 e 32- 2r? =0 r? =16 x= +4 


El valor x = —4 no interesa ya que está fuera del intervalo admisible. Por tanto, el único punto crítico de f en 
(0, +00) es z = 4. 
Puesto que 1? > 0 en (0, +00), es evidente que 


f'(x) > 0 en (0, 4) EN f es creciente en (0, 4) 
f'(x) <0 en (4, +00) f es decreciente en (4, +00) 


y que, en consecuencia, f tiene un máximo local en x = 4, que también es máximo absoluto de f en (0, +00). 


Solución: la distancia entre árboles óptima es de 4 metros. 


3. Se necesita encargar la fabricación de unas cajas de cartón para transportar plantas. Deben tener la forma 
que se indica en la figura, siendo la base un triángulo equilátero, y un volumen de 6 litros (6000 cm*). Se 
desea minimizar la cantidad de cartón necesario para la fabricación, esto es, el área total de la superficie de 
la caja, incluyendo las dos tapas. 


a) Deducir la función que proporciona el área total de las cajas que cumplen la condición dada (volumen 
igual a 6000 cm?) en función, sólo, de la variable z. ACLARACIÓN: no es la función del apartado b). 


b) Calcular la solución del problema 


Minimizar A(x) = > + — 
T 
para x € (0, +00) 


a) Comenzamos por escribir las expresiones del volumen del recipiente y de su área superficial. 


Volumen = V (x,y) = área base x altura 
Área = A(x, y) = 2 x área base +3x área lados rectangulares 


La base de la caja es un triángulo equilátero de lado x. Para cal- 
cular su área necesitamos conocer su altura, que se deduce usando 
el Teorema de Pitágoras: 
h 
2 3x? VB zx 


£ 
4 4 2 


2 ay? 2 2 2 


$ $ x/2 
En consecuencia se tiene 
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b) 


V3x V3a? me y V3 1*y 
V3 q? V3 2? 


Área rectángulo = z - y A(x,y) =2- J + 3xy = 3 + 3xy 


El volumen de las cajas debe cumplir la condición V (x,y) = 6000 cm y esta condición impone una 
relación de dependencia entre las variables del problema; por ejemplo, se puede despejar la y en función 
_4x6000 24000 


de la z: Ba 
37 
2=0000 + y= = 
4 V3 r? V3 r2 
y sustituyendo esta expresión de y en función de x en la expresión del área superficial tendremos ésta 
expresada, sólo, en función de la variable x, como nos piden: 


E 1 
Área triángulo = 21 


V3 q? V3 q? 24000  vV3x? 2400043 v3/ ə 48000 
A= + 31y = H3x. = = 6 ) 
2 2 BE 2 T 2 T 
jos 0 r? 24 ] 
Se trata de calcular el mínimo absoluto de la función A(x) = 5 + — en el intervalo (0, +00). 
zx 


Al tratarse de un intervalo abierto (que no contiene a sus extremos) ese mínimo, de existir, sólo puede ser 
un mínimo relativo (local). 


Estudiemos el crecimiento /decrecimiento de la función A. Su derivada primera es: 


2 24.24 2-2 
2 


T 


que sólo se anula para z = /24 = 2V3. 

Estudiemos ahora el signo de A’ en los intervalos (0,243) y (2V3, +00): 

En (0,243) se tiene z < 2/3 => z? < 24 > z? — 24 < 0, de donde 4” es negativa en (0, 2v3). 

En (2V3, +00) se tiene z < 243 => z3 > 24 > r? — 24 > 0, de donde A’ es positiva en (243, +00). 


Se tiene, pues, 
A'(x) < 0 en (0, 2/3) A es decreciente en (0, 2/3) 
=> 
A'(x) > 0 en (2v3, +00) A es creciente en (2V3, +00) 


= A tiene un mínimo local en z = 2V3, 


que, debido a la monotonía de A en (0,24/3) y en (2/3, +00), es, además, mínimo absoluto en el intervalo 
(0, +00). 


4. El dueño de una piscifactoría ha determinado que si compra x peces (en millares), entonces, al cabo de un 
mes tendrá 


a 9x 
—20+4 


f(x) 


peces. ¿Qué número de peces debe comprar para conseguir que la ganancia, f(1) — x, sea máxima? 


La cantidad óptima de peces que debe comprar será la que produzca el mayor valor de f(x) — x, es decir, la 
que produzca el máximo absoluto de la función f(x) — x en el intervalo [0, +00), ya que el número de peces a 
comprar será no negativo (en principio, cabe la posibilidad de que la máxima ganancia se produzca no comprando 
ninguno...) 

Luego el problema a resolver es 


Maximizar g(x) = f(x)— x 
para x € [0, +00) 
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g es continua y derivable en (0, +00). 
Derivada: 
(5 — 4x)(2x + 4) — (5x — 2z?)2  —4r? — 16r +20 _ 


g (2) = Qx +4) T rte 


Puntos críticos de g: 
g'(£)=0 > —4r? — 16r + 20 = —4(x? +4r — 5)=0 > r? +4r-5=0 


—4 + y 42 + 20 —4 +6 1 —5 


El valor x = —5 no interesa ya que está fuera del intervalo admisible. Por tanto, el único punto crítico de f en 
(0, +00) es z = 1. 


Conocidas sus raices, g'(x) se puede escribir g'(x) = —4(x + 5)(x — 1) y puesto que (2x + 4)? > 0 en (0, +00), 
es evidente que 


g'(x) > 0 en (0,1) g es creciente en (0, 1) 
g'(x) <0 en (1, +00) y es decreciente en (1, +00) 


y que, en consecuencia, g tiene un máximo local en x = 1, que también es máximo absoluto de g en [0, +00). 


Solución: el número óptimo de peces es 1 millar. 


5. La población de una especie sigue la siguiente función 
t 
P(t)=a+ 5, t>0 


Ea 


donde P(t) es el número de individuos de la población (medida en miles) y t el tiempo (medido en meses) y 
a es una constante positiva. 


) Calcular a sabiendo que inicialmente había 3000 individuos. 

) ¿En qué momento alcanza la población un máximo? ¿Cuánto es el valor de dicho máximo? 
c) ¿A qué tiende la población en el futuro? 

) 


Si se sabe que una población está en peligro de extinción cuando el número de individuos es menor que 
1000, ¿tiene esta población peligro de extinción? 


a) Tenemos que ver para qué valor de a se tiene 


P(0) = 3 (3000 individuos) =>  Ja=3 


b) El máximo absoluto de P(t) = 3 + a 3 + te™*/2 en el intervalo [0, +00) sólo puede ser un máximo 
e 


relativo o el punto t = 0. 


Veamos si P tiene algún máximo relativo: 
tn o 2 l e baja z 
Claramente se tiene, puesto que e7*/2 > 0 Yt € R, que 


P'(t) > 0 en [0,2) = P es creciente en [0, 2) 
P'(t) < 0 en (2, +00) => P es decreciente en (2, +00) 


Luego P tiene un máximo relativo en t = 2 que, claramente, es también máximo absoluto en [0, +00). 


El máximo absoluto de P en [0, +00) se alcanza en t = 2 y P(2) = 3.736 (3736 individuos) 


=h 
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c) Para ver a qué tiende la población tenemos que calcular 


t 
ím P(t)= lím 3+ 


t= +00 t= +00 et/2 


qu =3 


lm —>5 
t>+00 et/2 


lo que significa que | la población tiende a estabilizarse en 3000 individuos |. 


d) Obviamente, no hay peligro de extinción: 


> P(0) = 3 y P es creciente entre t = 0 y t = 2. 
> P es decreciente en (2, +00), pero no desciende del valor 3, al que tiende asintóticamente. 


Es decir, la población no desciende de 3000 individuos. 


6. La población de una especie sigue la siguiente función: 


(t =a)? 
Pa 


P(t)=1+ 


donde P(t) es el número de individuos de la población (medida en miles) y t el tiempo (medido en meses). 


(a) Calcula a sabiendo que inicialmente había 1700 individuos, esto es P(0) = 1.7. 
(b) ¿En qué momento la población aumenta? ¿Cuándo disminuye? ¿Cuándo alcanza un mínimo? 


(c) ¿A qué tiende la población en el futuro? 


a) Tenemos que ver para qué valor de a se tiene 


P(0) =1.7 (1700 individuos) 


17=1+3 a = a => 17+17=14+2%2 => 07=030 => |a= e 1.53 
b) Para ver cuándo la población aumenta/disminuye tenemos que estudiar el signo de P*(t): 
Pt) = 2(t— a)(1 + (t— ay?) RON _ Ata) +2 ay? 1 ay? __ At-a) : 
(1 + (t- ay?) (1 + (t — ay?) (1 + (t — a)?) 


Se ve claramente que 


e <0 para0<t<a > e decrece en [0, a) 


P'(t) >0 para t >a P crece en (a, +00) 


Obviamente, P tiene en t = a un mínimo relativo que es también su mínimo absoluto en [0, +00). 


—b= 
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c) Tenemos que calcular 


ta? 2 
E ad oeiia 


t— 2 
lím P(t)= lím ato ím 
t>o+o ] + (t — a)? t>+o0 12 


1 =14 
t>+o0 t>+o0 + 1+ (t =- a)? 


Esto significa que la población tiene a estabilizarse en 2000 individuos 


7. La población de cierta especie de pez que habita en un lago está dada por la siguiente función: 


i a 
POSTAT e t>0, 


donde P(t) es el número de individuos de la población (medido en miles) y t es el tiempo (medido en meses). 
Se pide: 
a) Calcula el valor de A sabiendo que inicialmente el lago tiene 2000 peces. 


b) ¿En qué momento la población de peces del lago alcanza su valor máximo? ¿Cuál es el valor de dicho 
máximo? 


c) ¿A qué tiende la población en el futuro? ¿En algún momento la población desciende por debajo de la 
población inicial? 


d) Esboza la gráfica de la función. 


a) Tenemos que calcular para qué valor de A se tiene 


P(0) = 2 (2000 individuos) 4> [A=2 


b) El máximo absoluto de P en [0, +00) sólo puede producirse en un máximo relativo o en el punto t = 0. 
Veamos si P tiene algún máximo relativo: 


Lia 
e == (1-,)=0 <= t=10 


1 


1 
Puesto que 9 e-w >0 para todo t € R, se tiene claramente 


t = 10 es un máximo relativo de P. 


P'(t) > 0 sit<10 
P'(t) <0 sit > 10 


Pero también es su máximo absoluto en [0, +00), ya que P es creciente en [0, 10) y decreciente en (10, +00). 


El máximo absoluto de P se obtiene en t = 10 y vale P(10) = 3.11 (3110 peces) 
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c) Tenemos que calcular 


o P t ¡O , t A 7 
lLím P(t)= lím (2+5e )=2+ lm ze s =2+ lm g gm0" 


t>+00 —>+00 


La función P(t) tiende asintóticamente a 2, lo que significa que la población tiende a estabilizarse en 2000 
peces y no desciende de este valor. 


d) 


8. La población de cierta especie que vive en un hábitat protegido sigue la siguiente función: 


t+1 
12424” Fa 


P(t) =A+2 


donde P(t) es el número de individuos de la población (medido en miles) y t es el tiempo (medido en años). 
Se pide: 

a) Calcula el valor de A sabiendo que inicialmente hay 3000 individuos. 

b) ¿En qué momento la población alcanza su valor máximo? ¿Cuál es el valor de dicho máximo? 


c) ¿A qué tiende la población en el futuro? ¿En algún momento la población desciende por debajo de la 
población inicial? 


d) Esboza la gráfica de la función. 


a) Utilizamos que P(0) = 3 (3000 individuos) para calcular A: 


2 2 70 
P(0)=A+5=3 = A=3- 57 <> A= 77 92.92 


b) El máximo (absoluto) de P(t) en el intervalo (0, +00), que es abierto, sólo se puede producir en un máximo 
relativo. Veamos si P(t) tiene alguno: 

t=-6 

t=4 


Pa) = P + 24— (t+1)2£ _ t? + 2t— 24 
(t? +24)? (t? + 24)? 


=0 => 4+2-24=0 => 


Vemos que el único punto crítico en el intervalo (0,00) es t = 4. Veamos si se trata de un máximo relativo 
estudiando el signo de la derivada a la izquierda y a la derecha del punto t = 4: 

t? + 2t — 24 
El signo de P'(t) = —2 ASNO es el opuesto al signo de t? + 2t — 24 = (t + 6)(t — 4), cuya gráfica es 
una parábola convexa. 
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t^2+2*t-24 —— 


20 e — 


Por consiguiente, 
t? + 2t— 24 < 0 en (0, 4) E P'(t) > 0 en (0, 4) z P es creciente en (0, 4) 
t? + 2t— 24 > 0 en (4, +00) P'(t) < 0 en (4, +00) P es decreciente en (4, +00) 


lo que demuestra que P tiene un máximo relativo en t = 4. Este máximo relativo es también máximo 
absoluto debido a que P es creciente en (0, 4) y decreciente en (4, +00). 


Luego, efectivamente, | P tiene un máximo absoluto en el punto t = 4 


El valor de P en dicho punto es | P(4) = > x 3.124 | 


c) Para saber a qué tiende la población en el futuro calculamos 


f t+1 E O 
jm (442777) =4+2 [im 777 =4=292 


Es decir, | La población tiende asintóticamente 2.92 (2920 individuos) 


La población inicial es de 3000 individuos, es decir, P(0) = 3. Hemos visto que, cuando t > 00, P tiende 
asintóticamente a 2.92. Lo cual significa que que, para t suficientemente grande, P toma valores tan 
próximos a 2.92 como se quiera. Obviamente, hay puntos en los que P toma valores menores que 3. 


La población desciende del valor inicial |. 


d) 


8 28 48 68 88 1089 128 148 


9. La población de una especie de insectos, tras la eclosión de los huevos puestos, sigue la siguiente función 


PO =r re, 
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donde t > 0 (t = 0 es el momento de la eclosión), siendo P(t) el número de individuos de la población 
(medida en miles) y t el tiempo (medido en días). 


a) 
b) 
c) 
d) 


Calcula a sabiendo que para t = 1 había 700 individuos. 
¿En qué momento alcanza la población un máximo? ¿Cuánto es el valor de dicho máximo? 
¿A qué tiende la población en el futuro? 


Esboza la gráfica de la función. 


a) Utilizamos que P(1) = 0.7 (700 individuos) para calcular a: 


b) 


e) 


d) 


P()=e"=0.7 —a = ln(0.7) = |a = —1n(0.7) ~ 0.36 


El máximo (absoluto) de P(t) en el intervalo [0, +00) sólo se puede producir en un máximo relativo o en 
t = 0 (extremo del intervalo incluido en el mismo). Vamos a calcular los máximos relativos (si existen) de 
P(t): 
t=0 
P' (t) =2te — ate = (2 — at)te7™™ =0 = 49 bien 3 


2-at=0 => t=7 5500 


Estudiamos el signo de P*(t) para saber dónde P es creciente/decreciente: puesto que e“ es siempre 
positivo, el signo de P'(t) coincide con el signo de (2 — at)t = 2t — at?, cuya gráfica, puesto que a > 0, es 
una parábola cóncava que se anula en t = 0 y en t = 2/a = 5.6: 


Por consiguiente: 


P'(t) > 0 en (0,5.6) 2 P es creciente en (0, 5.6) 
P'(t) < 0 en (5.6, +00) P es decreciente en (5.6, +00) 


En consecuencia, P tiene un máximo relativo en t = 5.6. 


P(5.6) = (5.6)? e7™36x5-6 ~ 4.18 > 0 = P(0) 


Está claro, pues, que | t = 5.6 es el máximo absoluto de P en [0, +00) |y que| P(5.6) = 4.18 |. 


Para saber a qué tiende la población en el futuro calculamos 


lím t e7% = (00 x 0) = lím z = (LH) =0 


t>00 t>o e 


Es decir, P(t) tiende asintóticamente a cero cuenado t > +00, lo que significa que 


la población tiende a la extinción |. 


— 10 = 
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10. Representar gráficamente la función siguiente, estudiando previamente dominio de definición, cortes con 
los ejes, asíntotas, crecimiento / decrecimiento, máximos y mínimos, concavidad / convexidad y puntos de 
inflexión: 


Dominio de definición: D = (—-oo, —1) U (—1, +00) 

Corte con el eje OY: f(0) = 1, luego la gráfica corta al eje OY en (0, 1). 

Corte con el eje OX: e” es siempre positivo, luego la gráfica no corta al eje OX 
Fu) <0six+1<0 es decir, si x < —1 

Fa) >0sizx+1>0es decir, si z > —1 


Signo de la función: f 


Asíntotas horizontales: y a 


Es decir, y = 0 es asíntota horizontal de f, para z —> —00. 
Asíntotas verticales: la única posible asíntota vertical es la recta x = —1, ya que este es el único punto de 
discontinuidad de la función: 


e” , e 


= 1m = 
+1 a>( + +1 


T 


lím 
a>(-1)- 


+00 


y) 


Asíntotas oblicuas: sólo podría haberla para x > +00, ya que para z > —o0 hay una asíntota horizontal. 
Si existe, será una recta y = mx +n con 


T 
m= lím Ft) = ím —&— = +œ = no hay asíntota oblicua. 
r=>+o g£ xr—> +00 glz + 1) 
Derivada: z 
iaon o te 
f (=)= (a + 1) 


que solo se anula para x = 0 que es, por lo tanto, el único punto crítico. 
Crecimiento y decrecimiento: estudiamos el signo de la derivada. El denominador de f'(x) es siempre 
positivo, luego el signo de la derivada coincide con el signo de xe”, que, a su vez, coincide con el signo de z: 


f'(x) <0 en (00, —1) f es decreciente en (—oo, —1) 
f'(x) < 0 en 1,0) y por lo tanto que f es decreciente en (—1, 0) 
f'(x) > 0 en (0, +00) f es creciente en (0, +00) 


=j 
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Extremos: como consecuencia de lo anterior, puesto que f es continua en x = 0 y pasa de ser decreciente a ser 
creciente en dicho punto, se tiene que f tiene un mínimo en xv = 0, en el cual f(0) = 1. 
Derivada segunda: 
Fa) (ze? +e%) (x + 1)? — 2xe* (2 +1) o (£ +1) — 2r e7 xr? +1 
Tt) = = e = 
(x +1)* (x + 1)? (x +1) 


que no se anula en ningún punto, ya que tanto e” como x? + 1 son siempre positivos. 


Convexidad y concavidad: estudiamos el signo de la derivada segun- 
da. Puesto que el numerador e” (g? + 1) es siempre positivo, el signo de a 
F"(x) coincide con el signo del denominador, (x + 1)?. Así, se tiene 


f” (x) <0 en (—00, —1) f es cóncava en (—oo, —1) 
F"(x) > 0 en (—1, +00) E f es convexa en (—1, +00) 


No hay puntos de inflexión. 5 > 


11. Representar gráficamente la función siguiente, estudiando previamente su dominio de definición, los cortes 
con los ejes, las asíntotas, los intervalos de crecimiento y decrecimiento, máximos y mínimos, intervalos de 
concavidad y convexidad y puntos de inflexión: 


l1+lnz 
E 


f(x) 


T 


Dominio de definición: D = R+ = (0, +00) 

Corte con el eje OY: no hay, ya que la función no está definida en x = 0. 
Corte con el eje OX: f(1)=0 $ l1+Inz=0 $ r=e!. 

f(z) <0si1+lng < 0 es decir, si z < e7! 

f(x) > 0si 1+Inx > 0 es decir, si x > e”? 


Signo de la función: i 
Asíntotas horizontales: 
1+Inx [LH] 1/x 1 


ím —T—= lím =>=0 
2 z>+00 2g r=>+o 212 


lím 
zT=>+00 XT 
Es decir, y = 0 es asíntota horizontal de f, para x —> +00. 
Asíntotas verticales: la única posible asíntota vertical es la recta x = 0 
l+lnz  —oo 


lím = = 00 
x>0+ q? 0 


Asíntotas oblicuas: no hay, ya que para x —> +00 f tiene una asíntota horizontal. 
Derivada: iasi 
r —-1-2Inzx 

2) = ——— 
re == 
que solo se anula para x =e71/2 que es, por lo tanto, el único punto crítico. 
Crecimiento y decrecimiento: estudiamos el signo de la derivada. El denominador de f'(x) es siempre 


positivo en (0, +00), luego el signo de la derivada coincide con el signo de —1 — 21n z: 


F'(x) > 0 en (0,e7?/?) 
F'(x) <0 en (e712, +00) 


f es creciente en (0,e71/?) 


lo t 
DOMO ds i f es decreciente en (e71/?, +00) 


Matemáticas Aplicadas a la Biología Colección de ejercicios resueltos. Tema 1 


Extremos: como consecuencia de lo anterior, puesto que f es continua en x = e71? y pasa de ser creciente a 
ser decreciente en dicho punto, se tiene que f tiene un máximo en z = e72, 
Derivada segunda: 


2 
3 2 
2% + (1 + 2lnz)3z -2+3(142Inx) 1+6lng 


1 = T == 
f (x) a xê pl xt 4 


1/6. 


z 
que sólo se anula si 1 + 6 ln x = 0, es decir, para z = e7 
Convexidad y concavidad: estudiamos el signo de la derivada segun- „4 


da. Puesto que el denominador x4 es siempre positivo, el signo de f”(x) 
coincide con el signo de su numerador, 1 + 6ln x. Así, se tiene 


f'(x) < 0 en (0,e7 1/8) 3 f es cóncava en (0, e76) 
f"(x) > 0 en (e716, +00) f es convexa en (e, +00) Le 


= 3 X 
ee ¿16 


Obviamente, en x = e716 la función tiene un punto de inflexión. 


12. Representar gráficamente la función siguiente, estudiando previamente su dominio de definición, los cortes 
con los ejes, las asíntotas, los intervalos de crecimiento y decrecimiento, máximos y mínimos, intervalos de 
concavidad y convexidad y puntos de inflexión: 


SS 
AOS 


Dominio de definición: D = R, ya que e”! no se anula. 


Corte con el eje OY: f(0) = z5 3e, luego la gráfica corta al eje OY en el punto (0, 3e). 
Corte con el eje OX: f(1)=0 & r?-—4r+3=0 & r=3ór=l1. 
Signo de la función: el denominador es siempre positivo y el numerador puede escribir factorizado 
f(x) > 0 en (oo, 1) 
x? — 4x + 3 = (x — 3)(x — 1), de donde se deduce facilmente que 4 f(x) < 0 en (1,3) 
f(x) > 0 en (3, +00) 


Asíntotas horizontales: 
, z?—4r+3 +œ 2  a2—4di+3 +o 
lím = =+00 y lím = 
2>—00 ex-1 0 u>+o0 ex-1 +00 


= (LH 2 veces) = 0 


Es decir, f tiene una asíntota horizontal en y = 0, para x —> +oo. 
Asíntotas verticales: no hay, ya que f es continua en R. 
Asíntotas oblicuas: no hay. 
Derivada: 
Fo (2x — 4) eel (x? —4x4+3).er1 (2g — 4) — (x? — 4x + 3) r? +6r-—7 
t) = = = = 


(er-1y? er—1 er—1 


que se anula si —x? + 6x — 7 = 0, esto es, para x = 3 — v2 y para z = 3 + v2. 
Crecimiento y decrecimiento: estudiamos el signo de la derivada. El denominador de f'(x) es siempre 
positivo. El polinomio del numerador cambia de signo en x = 3 — y2 y en z = 3 + y2, y se puede ver que 


F' (a) < 0 en (—00,3-— v2) f es decreciente en (—00,3— v2) 
f'(x) > 0 en (3— v2,3+ v2) y por lo tanto que f es creciente en (3 — V2,3 + v2) 
F'(x) < 0 en (3 + v2, +00) f es decreciente en (3 + v2, +00) 


TB 
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Extremos: como consecuencia de lo anterior, se tiene que f tiene un mínimo local en z = 3 — y2 y un máximo 
local en z = 3 + v2. 
Derivada segunda: 


z? — 8z +13 
f” (£) = RF 


que se anula en z = 4 — v3 y en z = 4 + v3. 


Convexidad y concavidad: estudiamos el signo de la de- 
rivada segunda. Puesto que el denominador es siempre po- 
sitivo, el signo de f” (x) depende del signo de su numerador, 
x? — 8x + 13, y se tiene 


f"(x) > 0 en (=, 4 — v3) 
f'(x) <0 en (4— V3, 4 + v3) 
f(x) >0en (,4+ v3, +00) 


f es convexa en (—00,4— v3) 
f es cóncava en (4— v3, 4 + v3) 
f es convexa en (,4 + v3, +00) 


Obviamente, en z = 4 — V3 y x£ = 4 + v3 son puntos de 
inflexión. 


Ejercicios del Tema 2: Modelos discretos 
en Biología 


13. La dinámica de una determinada población de mamíferos (medida en miles) que habita un parque natural 
viene dada por el modelo de Ricker: 


dba = mad e), n>1 


con parámetro de crecimiento r = 0.4 y capacidad de alojamiento K = 5 (miles). Se pide: 


a) Si la población inicial de mamíferos es zo = 3.5, ¿cuál es la población tras dos periodos de reproducción? 
b) Calcula los puntos de equilibrio del modelo y estudia su estabilidad. 


c) Interpreta biológicamente los resultados. 


1. Calculamos z1 y t2, a partir de £o = 3.5: 
21 = roe UR) = 3.50%) e 3.946 
zr = rye" UR) =3.9460010-*89) y 4,293 
Esto es, tras un periodo reproductivo habrá 3946 ejemplares y después de dos habrá 4293. 


2. Los puntos de equilibrio son las soluciones de la ecuación z = f(x) con f(x) = xer027/K) 


z=0 
a = qe (0-7/K) == z(1— er(oz/K)) =0 => o bien 
erli2/0) 1 => 1-¿=0 => r=K 


Los puntos de equilibrio son ri = 0 y z$ = K = 5 


Para estudiar la estabilidad de los puntos de equilibrio calculamos en primer lugar la derivada de f 


f(x) = er(l-2/K) y a Ja 2 a T EAR 


> f(x) = f'(0)= (1-0)e” = e” > 1 (porque r > 0), luego el punto de equilibrio| xj = 0 es inestable 


2 
> f (23) = (5) = (1- El er1-5/5) = a luego el punto de equilibrio | 13 = 5 es estable 


3. D xf = 0 es un punto de equilibrio inestable: Esto significa que si inicialmente hay 0 individuos, seguirá 
habiendo 0 individuos para todo tiempo posterior. El hecho de que sea inestable indica que, si se 
inicia el modelo con un número zo Æ 0, aunque sea próximo, el sistema evolucionará alejándose de 0. 


15 
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> 23 = 5 es un punto de equilibrio estable: Esto significa que, si inicialmente hay 5000 individuos, la 
población se mantiene en equilibrio en dicho valor (recuérdese que esa es la capacidad de alojamiento). 
Este punto de equilibrio es estable, lo que significa que si se comienza con un valor próximo a K = 5, 
aunque distinto, el sistema evolucionará acercándose al punto de equilibrio. 


14. La dinámica de una población de peces sometidos a una pesca severa sigue el siguiente modelo discreto 


Epa = Bo + 2% (1 — =; — 1, 


donde x,, es el número de peces (medido en miles) en el año n. 


a) Si zo = 2, ¿cuál es el valor de x2? 
b) Calcula sus puntos de equilibrio y estudia su estabilidad. 


c) Interpreta biológicamente los resultados. 


2 
. 21 =%0+2m(1- F) -1=24+4(1-7)-1=3 
3 
z =a +20 (1-7) -1=3+6(1-3)-1=35 


2. Tenemos que resolver la ecuación x= f(x) con f(x) = x +2x(1— 3) 1: 


2=2-— V2 x 0.58 
— 4z +2 = 0 < o bien 
r =2 + 2 x 3.41 


sr=r+2æ(1-)-1 < 20(1-5)-1=0 => q 


Los puntos de equilibrio son xj = 0.58 y 13 = 3.41 (valores aproximados) |. 


Para estudiar la estabilidad de los puntos de equilibrio necesitamos calcular la derivada de f: 
Flo) =3-x. 


Se tiene: 


F (21) = 3 — 0.58 = 2.48 > 1 => 21 = 0.58 es localmente inestable 


F (23) = 3 — 3.41 = —0.41; |f (23) < 1 => | 13 = 3.41 es localmente estable 


. Desde el punto de vista biológico, que zł = 0.58 y 13 = 3.41 sean puntos de equilibrio significa que si el 
sistema se inicia con estos valores, permanece constante, es decir el número de individuos de la población 
no cambiará en el tiempo. En efecto, si se comienza con zo = 0.58 (miles de peces), 


0.58 
xo = 0.58, 11 = £o + 2xo(1 — D) — 1 = 0.58 +2 x 0.58(1 — Po — 1 = 0.58, etc. 


Análogamente, si se comienza con £o = 3.41 (miles de peces), 


3.41 
1 =3.41, 2 =20+20(1-7)-1=3414+2x341(1- =) -1=341, etc. 


Que el punto de equilibrio x = 3.41 sea localmente estable quiere decir que, si iniciamos el sistema con 
un valor xy Æ 3.41, pero próximo, el sistema evolucionará acercándose a este valor. 


Por el contrario, que el punto de equilibrio 13 = 0.58 sea localmente inestable significa que, si comenzamos 
con un valor xy Æ 0.58, por muy próximo que sea, el sistema evolucionará alejándose del valor 0.58. 


— 16 = 
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15. El tamaño de una población de insectos (medida en millones) sigue una dinámica dada por la curva de 


reclutamiento de Beverton-Holt: 
il 
my = == 


A 
K Tk-1 


con constante de crecimiento R = 1.4 y capacidad de alojamiento K = 10 (millones). Se pide: 


a) Si la población inicial de insectos es £o = 4.5, ¿cuál es la población tras dos periodos de reproducción? 
b) Calcula los puntos de equilibrio del modelo y estudia su estabilidad. 


c) Interpreta biológicamente los resultados. 


1. Pongamos 


Ko Rz 1.4x l.4x 
LD) = = = 
R-1 1.4-1 1+0.04x 
1 
A A 
A partir de xy = 4.5, calculamos zı y La: 
las e 2800 
tE = = k = ————— 7x . 
E= AMO 1 +0.04 x 4.5 
1.4 x 5.339 
= = f(5.3390) = —————————— 7x 6.1592 
z2 = f(21) =$ ) = 140.04 x 5.339 
2. Tenemos que resolver la ecuació f(x) con f(x) sa 
. Tenemos que resolver ión £ = n = 
q i A 7 = 170.042 
g = 
1.4x 1.4 ; 
| 4 1%) =0 = q o bien 
“140.047 o 1+0.042 14, sn 
1+ 0.047 dl 


Los puntos de equilibrio son xi = 0 y z$ = 10 | 


Para estudiar la estabilidad de los puntos de equilibrio calculamos la derivada de f: 


Ps) = 1.4(1 + 0.04 7x) — 1.4 x 0.04% 1.4 


(1 + 0.04 2)? (1 + 0.04 2)? 


Se tiene: 


Far) = f'(0) = 1.4 > 1 => | xf = 0 es inestable 


1 
F (23) = (10) = I~ 1 => 23 = 10 es estable 


3. La interpretación biológica de estos resultados es la siguiente: 


Si se inicia el sistema con xy = 11 = 0 individuos, obviamente se mantiene constante igual a cero: no hay 
población. Pero, que sea un equilibrio inestable significa que si comenzamos con un valor de xq distinto 
de 0, aunque sea próximo, el sistema evolucionará alejándose del valor 0. 


Por otra parte, si iniciamos con zo = 13 = 10 (millones de individuos) también se mantendrá constante e 
igual a esa cantidad. Pero, como además es estable, si comenzamos con otro valor de zo, el sistema tenderá 
asintóticamente a ese equilibrio, es decir, la población tenderá asintóticamente a 10 millones. 


== 
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16. En una cierta colonia de focas, cuyas hembras se han clasificado en 3 grupos de edad, se han observado las 
siguientes tasas de supervivencia: un 60% de las hembras de edad 0 sobrevive hasta la edad 1, mientras que 
un 85% de estas últimas lo hace hasta la edad 2. Ninguna sobrevive más allá de la edad 2. 


Por otra parte, se ha estimado que las hembras de edad 1 tienen una tasa de fecundidad del 80% (es decir, 
tienen de media 80 cachorros por cada 100 hembras), mientras que las de edad 2 tienen una tasa del 50%. 
a) Escribir la matriz de Leslie asociada a estos datos. 


b) Partiendo de una población inicial de 100 focas de edad 0, de 60 de edad 1 y de 50 de edad 2, calcula 
el número de focas de edad 2 que habrá tras 2 períodos de observación. 


c) En una visita a la colonia se recuentan 85, 85 y 60 focas de edades respectivas 0, 1 y 2. ¿Qué distribución 
por edades se puede suponer que hubo en el período anterior? 


1. Denotemos E?, E! y E? a las tres clases de edad. Los datos que nos dan se reflejan en el siguiente esquema: 


t — a] E 


Escribamos estas relaciones en forma de ecuaciones y en forma matricial: 


Ep, =0.8E; + 0.5Ef E? 0 0.8 0.5 E? 
Era = 0.6£? EL, |=||06 0 0 El 
Efy1 =0.85E; Eres 0 0.85 0 E? 


2. Calculemos, según el modelo anterior, la distribución de la población en los instantes de tiempo k = 1 y 


a E? 0 08 0.5\ /El 0 0.8 0.5\ /100 73 
El =ļ|o6 o oļļzļ|=ļo0o6 o o|] 6so]ļ=]l60 
E? o 0.85 0) 2 0 0.85 0 50 51 
El 0 08 05 [El 0 08 0.57 /73 73.5 
El|=|o06 o 0 ]||eÉi]|=lfo6 o 0 ||60|=|438 
E2 o 0.85 0/ 2? o 0.85 0/ \51 51 


Tras 2 periodos de observación habrá 51 focas de edad 2 


3. Supongamos que la visita a la colonia se produce en el instante de tiempo k, es decir, que 


E? 85 
El |= |85 
E? 60 


Siguiendo el modelo, se tiene la siguiente relación entre la distribución de la población en el instante k y 
en el instante k — 1: 


=j 
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E? 85 0 08 05\ (E0, 
El|=|85s|=|o06 0 0 ||E£l, 
E? 60 o 0.85 0/ AE, 


Es decir, la distribución de edades del periodo anterior viene dada por la solución de este sistema lineal 
de ecuaciones. De la última ecuación se tiene 


60 6000 


1 1 
0.85Ej_1 =60 —= Epi = a = 70.582 
De la segundo ecuación: 
85 850 
06, ¿== El = 0676 1416667 


Y de la primera ecuación: 


0.8E%_1 +0.5E%_, =85 <= 0.5E%_, =85-—0.8E;_, = 85 — 0.8 x 70.5882 = 28.5294 
28.5294 


= 57.0588 
0.5 


- = 
<> B= 


Así pues, en el periodo de tiempo anterior al recuento, según el modelo, la población estaba distribuída 
de la siguiente forma 


E? 141.6667 
EL, | = | 70.5882 
El, 57.0588 


17. Supongamos que en una cierta población, cuyos individuos se han clasificado en 3 grupos de edad, la matriz 
de Leslie es de la forma 


2 AL: 
04 0 0 
0 0.7 0 


a) Interpretar cada elemento de esta matriz. 


b) Determinar lo que sucede si se realiza un seguimiento de una población durante dos estaciones repro- 
ductivas, comenzando con 1000 individuos de edad 0, con 100 de edad 1 y con 50 de edad 2. 


1. Denotemos E%, E! y E? las tres clases de edad. 


> La primera fila de la matriz de Leslie da cuenta de la reproducción de la población: 
e Por cada individuo de la clase E% nacen, en promedio, 2 nuevos individuos en cada estación 
reproductiva. 
e Por cada individuo de la clase E! nacen, en promedio, 5 nuevos individuos en cada estación 
reproductiva. 
e Por cada individuo de la clase E? nace, en promedio, 1 nuevo individuo en cada estación repro- 
ductiva. 


> La subdiagonal da cuenta de la supervivencia. En este caso indica que un 40% de los individuos de 
edad E? consigue sobrevivir hasta el periodo vital siguiente y que un 70% de los individuos de edad 
E! sobrevive hasta el siguiente periodo reproductivo. 


— 19 = 
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> El resto de elementos de la matriz, que son nulos, indican que no hay otras interrelaciones entre las 
clases de edad. 


2. Calculemos cuál será la distribución de los individuos en clases de edad en los dos periodos reproductivos 
siguientes, si partimos de la distribución inicial: 


El 1000 

E} | = | 100 

E? 50 
E? 2.5 1\ /E 2 5 1\ /1000 2550 
En|=|o4 o oļjļ|zi|=]ļo0o4 o o|| 100 |= f 400 
E? o 0.7 0) E o 07 0) A 50 70 
El 2. 5 mM/2 2 5 1\ /2550 7170 
El|=|o4 o o||E¡|=fo4 o o|| 400 | =| 1020 
E2 o 0.7 0/ \E? o 0.7 0/ \ 70 280 


18. Tenemos una población (de hembras) dividida en dos clases de edad, llamadas Et y E? respectivamente, y 


cuya matriz de Leslie es 
L= 2 di 
2 X.0.75 0)” 


(a) Interpreta biológicamente cada elemento de la matriz de Leslie. 

(b) 

(c) Discute el comportamiento en el futuro de la población y calcula las proporciones a largo plazo de las 
clases de edad. 


Calcula el autovalor dominante (principal) y un autovector asociado. 


1. El modelo de Leslie asociado a esta matriz es: 
ELN (2 A(E 
E2)  \075 0) |2} 


> La primera fila de la matriz de Leslie da cuenta de la reproducción de la población: 


200 
e Los individuos de la clase Æt se reproducen a una tasa de 2 = 100 ~ 200%, es decir, por cada 


100 individuos de la clase Æ! nacen, en promedio, 200 nuevos individuos en cada periodo de 
tiempo. 


Interpretación: 


400 
e Los individuos de la clase E? se reproducen a una tasa de 4 = 100 ~ 400 %, es decir, en promedio, 


por cada 100 individuos de la clase E? nacen 400. 


> La subdiagonal da cuenta de la supervivencia. En este caso indica que un 75% de los individuos de 
edad E! consigue sobrevivir hasta el periodo vital siguiente. 


> El 0 del último elemento de la diagonal significa que ningún individuo de la clase Æ? sobrevive más 
de un periodo vital. 


2. Los autovalores son las soluciones de la ecuación: 


= 2— À 4 S -49 a à =3 
det(L— A) = det (y; $) 0-M-3=1 po E 0 - 


=a jm 
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Vemos que à; = 3 es un autovalor dominante, ya que |3| > | — 1| = 1 y además es positivo. 


Calculamos ahora un autovector asociado al autovalor dominante A; = 3, es decir, calculamos una solución 
del sistema lineal de ecuaciones 


(L-31) (1) E (os 5) (a) p (o) 


Se trata (como ya sabemos, por ser À un autovalor) de un sistema compatible indeterminado: las dos 
ecuaciones del sistema son equivalentes. De una cualquiera de ellas se tiene: 


vı = da : 
—v + 4v = 0 ví = 4u» 1 para cualquier aœ € R 
V2 =Q 


Luego un autovector es, por ejemplo, (eligiendo «a = 1) | v = (i : 


„ Tkp 
3. Por teoría sabemos que lím a 
k>+00 Tk 


Esto significa que, asintóticamente, Tk+1 ~% 3Tk, es decir que, a largo plazo, la población se multiplica por 
3 en cada periodo de tiempo. En consecuencia, la población crece indefinidamente (porque 3 > 1). 


= A = 3 siendo Ty la población total en el instante k: Ty = E + Ef. 


Por otra parte, también sabemos por teoría que las proporciones de las clases de edad frente al total 


verifican 
El El 4 4 
A o “i =2=08 
Ty Ez + E% k=>00 V1 + V2 4+1 5 
E? E? 1 1 
Eoee a a E = 20.2 
Tp Ez + E; k=>00 UV1 + V2 4+1 5 


Es decir, a largo plazo, la clase de edad E! será el 80% de la población mientras que la clase de edad E? 
será el 20%. 


19. Tenemos una población dividida en dos clases de edad, llamadas E! y E? respectivamente. Se sabe que la 
cantidad de individuos nuevos en cada generación se obtiene a partir de un promedio de 0.2 veces el número 
de individuos de tipo E! y 0.6 veces el número de individuos de tipo Æ?. Por otro lado, en el paso de una 
generación a otra, la probabilidad de supervivencia de la etapa Æt hasta la etapa E? es del 25%. 


(a) Construye la matriz de Leslie asociada al modelo. 
(b) Calcula el autovalor (dominante) principal y un autovector asociado. 


(c) Estudia el comportamiento en el futuro de la población y calcula las proporciones a largo plazo de las 
clases de edad. 


1. La información sobre la reproducción y la supervivencia que nos dan se reflejan en el siguiente esquema: 


Escribamos estas relaciones en forma de ecuaciones y en forma matricial: 
a 0.2 0.6 El 
E2 J \o25 0J)Nel 


— lp 
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2. Los autovalores son las soluciones de la ecuación: 


0.2-A 0.6 
0.25  —A 


A1 = —0.3 


det(L — AI) = det ( 
à2= 0.5 


) = —\A(0.2 — A) — 0.15 = à? —0.24 — 0.15 = 0 => i 
Vemos que A2 = 0.5 es un autovalor dominante, ya que 0.5 > | — 0.3| = 0.3, y además es positivo. 


Calculamos ahora un autovector asociado al autovalor dominante A2 = 0.5, es decir, calculamos una 
solución del sistema lineal de ecuaciones 


(L—0.51) a = E e) A = (o) 


Se trata (como ya sabemos, por ser A3 un autovalor) de un sistema compatible indeterminado: las dos 
ecuaciones del sistema son equivalentes. De una cualquiera de ellas se tiene: 


=2 
—0.3v1 + 0.6v = 0 vi = 209 a para cualquier a € R 
V2 = Q 
; Eh 2 
Luego un autovector es, por ejemplo, (eligiendo a = 1) | v = 1)! 
Tr 
3. Por teoría sabemos que g lím - = M2 = 0.5 siendo Tķ la población total en el instante k: Te = E} + E2. 
—>+00 k 


Esto significa que, asintóticamente, Tķ+ı > 0.5 Tk, es decir que, a largo plazo, la población se multiplica 
por 0.5 en cada periodo de tiempo. En consecuencia, la población decrece a cero (porque 0.5 < 1), esto 
es, tiende a la extinción. 


Por otra parte, también sabemos por teoría que las proporciones de las clases de edad frente al total 


verifican 
El El 2 2 
Es E e =2 ~ 0.67 
Th Ex + Ez k=>00 vi + U2 2+1 3 
E? E? 1 1 
Bec a 2 = > x 0.33 
Ty Eg + Ex k=>00 Ur + U2 2+1 3 


Es decir, a largo plazo, la clase de edad E! será el 67% de la población mientras que la clase de edad E? 
será el 33%. 


20. Para estudiar la dinámica de una población de roedores, cuya vida media es de 6 meses, se divide la misma 
en tres grupos según su edad, E), EY) y EU), y se obtienen los siguientes datos sobre su fertilidad y 
supervivencia: Sólo un 17% de los individuos de edad E) llegan a la edad EY) y sólo un 24% de los de 
edad EU” llegan a la edad EG). Además, un 5% de los roedores de edad EP? consiguen sobrevivir un periodo 
extra. Por otra parte, en promedio, la fertilidad del grupo EU? es de 6 crías por individuo, mientras que la 
de la edad EG) es de 3. 


Escribe la matriz de Leslie asociada a estos datos. 


— 2 
id NS 


Un grupo formado por sólo 100 individuos de edad EG) ocupa un espacio en el que no había ningún 
roedor previamente. Calcula el número de individuos de cada edad que habrá en dicho espacio tras dos 
periodos de observación. 


c) En un recuento concreto, se observa que hay 900 roedores de edad EU), 34 de edad EG), y 29 de edad 
EGB). Calcula la distribución por edades que debió haber en el periodo anterior. 


d) Sabiendo que el autovalor principal, que es dominante, y un autovector principal de la matriz de Leslie 


son: 
0.99 


A»=107,  Us=|[0.16|, 
0.04 


— 0 — 
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calcula el comportamiento en el futuro de la población y la distribución estable de cada clase edad. 


0 6 3 
L=| 0.17 0 0 
0 0.24 0.05 


2. Los datos proporcionados corresponden a la condi- | Tiempo | Edad 1 | Edad 2 | Edad 3 | 
ción inicial EE zO EO | aO |] 
0 l f 
_ o 0 100 07 
di ri Ca 600 0 24 | 
E 72 102 | 12 | 


Calculamos P, y Pa aplicando dos veces la ecuación 
de recurrencia Py = LPk—1 


3. En este caso, en la relación de recurrencia entre Pk y Pk-1, Pk = LP;-1, conocemos Py y queremos 
calcular Pg. Para calcular el número de roedores en el período anterior, hay que resolver el sistema 


lineal: 
1 1 (2) (3) E 
Grea EN, EN i CE, +35 = 900 
017 0 0 |[|E8, |== |E® |= |34] = {0.172% =34 
O 0.24 0.05 29 
B®, El 0.ME%, +0.05E8%, =29 


De la segunda ecuación: 


De la primera ecuación: 


900 — 6E% 
GE, T So 900 E, 3 E=1 = 300 2 EO, 


De la tercera, sustituyendo ES: 


0.24 E® +0.05E%, =29 <= 0.24 E® +0.05(300— 2 E, ) = 29 
<= 0.14E?, +15=29 


<= E?, =100 


y, en consecuencia 


3 2 
EC =300-—2E(?, = 300 — 200 = 100 


Luego, la generación anterior se componía de (200, 100, 100} roedores respectivamente. 


— o 
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4. Como sabemos, por el Teorema 2.38, se tiene 


T 
ím —* 


=A 
k— oo Tk—1 p> 


lo que significa que, para valores grandes de k, se tiene Tk ~ Ap Tx-1; es decir, la población total, para 
tiempos grandes, se comporta según un modelo exponencial con constante de crecimiento R = Ap. Como 
Ap > 1, eso significa que la población total tiende a crecer indefinidamente. 


El Teorema 2.38 también nos proporciona información sobre el comportamiento asintótico de las fracciones 
de las clases de edad con respecto al total. Concretamente: 
E) 
lím 
k=>00 Tk 


B 0.99 
-0.99 +0.16 + 0.04 


= 0.8319, 


lo cual significa que, para tiempos suficientemente grandes, el grupo EU) será el 83.19% de la población 
total. 


Análogamente, para los grupos E) y E): 


ES) 0.16 
í z = 0.1345 
k>% T, 0.99 +0.16 + 0.04 ! 
EG) 0.04 
Y = = 0.0336 
k> T, 0.994 0.16 +0.04 ; 


lo cual significa que, para tiempos suficientemente grandes, los grupos E(% y EG) constituirán, respecti- 
vamente, el 13.45% y el 3.36% de la población total. 


21. Comprueba que (0,0) es un punto de equilibrio del modelo siguiente, y estudia su estabilidad: 


mp = YUk T 
Alark a) 
o 820 
= Tk—1 
IE 


Por definición, un punto de equilibrio es cualquier punto (x, y) que verifique: 


“== 1 
£ 
Y= 1 2 904) 


Sustituyendo en estas ecuaciones x = 0 e y = 0, resulta obvio que (0,0) es un punto de equilibrio. 


Para analizar su estabilidad, debemos calcular los autovalores de la matriz jacobiana J(0, 0). 


Comenzamos por calcular las derivadas parciales de las funciones f y g: 


ar ) ai ) —61y 1 

an NU) Y a At, Y p 

0) ð 9(1 +x?)  3(1+2? 

maot y _ [91 +=) (1+2?) luego J(0,0)= a an 
99 ) iy ) 1 —2xy 1 0 
de z, Yy Əy Y 1+ y (1 + y?) 


Calculamos ahora los autovalores de J(0, 0): 


det (J(0,0) — À Id) = det ( 


=A 
1 


Mi) =X- 1/3=0 <> A=zw 


— 
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Puesto que ambos autovalores son, en módulo, menores que 1, 


1 -1 
=—|=|=|< 1 


concluimos que el equilibrio nulo (0, 0) es estable. 


22. Calcula los equilibrios biológicamente relevantes del modelo Nicholson-Bailey y analiza su estabilidad. 


fs = Demo E RRA, 


Yk = 2 2-1 (1 — e lvn-1), 


Por definición, un punto de equilibrio es cualquier punto (x, y) que verifique: 


z = 20070 = f(z,y) 
y=2w(1—e7%Y) = g(x,y) 
Se trata de un sistema (no lineal) de dos ecuaciones con dos incógnitas. Vamos a resolverlo despejando alguna 


de las variables en una de las ecuaciones y sustituyendo en la otra. 
De la primera de las ecuaciones: 


zr=0 
a=200 Y => (1-2 =0 = o bien r , 
1= 2000 > 0000 y = —10 ln (5) = 10In(2) 


Vemos pues que la primera ecuación se verfica si x = 0, o bien si y = 101n(2). 


En el primer caso, de la segunda de las ecuaciones se tiene, sustituyendo x = 0: 


y = 2% (1 — poa =0. 


Es decir, | (x, y) = (0, 0) es un equilibrio |. Este equilibrio (trivial) no tiene ningún interés (no es relevante). 


Para explorar el otro caso, sustituimos y = 101n(2) en la segunda de las ecuaciones del sistema: 


y=20(1-e%'Y%) <= 10In(2) = 2x(1 — >) =% 


Luego, | (x,y) = (101n(2), 101n(2)) = (6.93, 6.93) es otro equilibrio |. 


Para analizar su estabilidad, debemos calcular los autovalores de la matriz jacobiana J(101n(2), 101n(2)). 


Ay) Loy) 
Ir Y Oy Y 2e70-ly —0.2xe 0-14 
Ja, y) a ð = —0.1 T 
Ir, y) L y) A 
dx” dy” 
luego, 1 
ge — ln(2) 
5) 1 —1In(2) 
J(101n(2), 10In(2)) = 1 = G mo) 


2(1-5) In(2) 


Calculamos ahora los autovalores de J(101n(2), 101n(2)): 


— 
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det (Q0 Mm(2), 1010(2)) — Ad) = det (17> a) 


= (1 — A)(ln(2) — A) + In(2) 
= M — (1+1In(2)A + 21n(2) = 0 


Las raíces de esta ecuación de segundo grado son: 


y — @tn(2)) + VQ +In(2)2-81In(2)  (1+I(2)) + y1 + ln? (2) — 4In(2) _ 1.69 + v=1.29 
o 2 o 2 2 
Como vemos, los autovalores son complejos conjugados. Utilizando la Observación 2.57 se tiene, puesto que 


|21n(2)| = 1.386 > 1, que ambos autovalores tienen módulo mayor que 1. 


Luego | el equilibrio no trivial (10 In(2), 10 In(2)) = (6.93, 6.93) es inestable |. 


— 6 


Ejercicios del Tema 3: Cálculo integral 


23. Calcular la siguiente integral indefinida 


1 
z? pan > sE ne)+ fa i 
z zx 
ln(x) i E 
> z7 C= (In(e) + 1) +C 


1 A B 


= 0 > A= 1 
1= A(z +1) + Bz: 4 
TA seai > (x + 1) + Bz; i 


1 1 1 
|z” JCG =) x = ln |z| — ln |x + 1| +C = ln 


luego, 


25. Calcula la siguiente integral indefinida: 


cos(t) sen(t) 
l (2 + cos(t))? a 


Esta integral se calcula haciendo un cambio de variable: 
J cos(t) sen(t 


OS la a ] = | arp” = Itri) EAS 


27 
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2 1 2 1 2 
ata Ju 22H u ENT o] 2 + cos(t) 
(*) Reducción a suma de fracciones simples: 
—u A B u=-2>2=B 
= | = A(2 B 
(+07 244 tu» e cs E ela 0=24+2>A=-1 
es decir, 
-u _ -l T: 2 
(2+u) 2+u (2+u) 
26. Calcular la siguiente integral: 
0 In(3x) 
3 dx 
de 
In(3x) u=1l0(31); u = : 1 1 
n(3x = ; == — 
= T = =—z l — = 
I z5 dx Je o9 „= 1 292 n(3x) +f 275 dx 
as” 22? 
o) = EN, aea In(3x) L pge a EU g 
2m2 TOT Ta a? ET T 4r? 
27. Calcular la siguiente integral: 
0 In(1 + 2?) de 
2x 
2 2 


u=In(1+232); u 


Jara?) de = | P 


28. Calcular la siguiente integral: 


2x 


—  de= 
| id 


1 


Jl 7) 


1; 


2 


oln(t+a)= f2- + 


F 
1422 


V=T 


| 


dx =xIn(1 +2?) — 2z + 2arctgr +C 


q? 
J e27 
1 + 3e7 
t= e 0) t a= f 1_ 1⁄3 
dt = edx | J 1+3 3 3t+1 
1 1 1 1 
dt== | t— - ln|l3t+ 1 =|= [e”--1 “+1 
¿ll 3 n |3t + |) +c 1C a n |3e? + +0 


—:98:= 
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(x) Obsérvese que e?” de = e” e” da = tdt. 


29. Calcular la siguiente integral: 


| lng d 
x(3+21nx) dl 


nz t=Inw t 1 3/2 
——— de = 1 = | — at= | [-- t= 
eo ij pen E IG “OL 
T 


1 3 1 1 1 3 
1 d=- |t In |2t y =|- (n(x) — ? 1/21 
Jl m3) 7 13 i 3l) +C 3 (mt) > n| na) +3) +0 


30. Calcular el área encerrada por las gráficas de las parábolas 


y= 2x — Tr +5 e y =-22 4827 


y = 2x? — Tx + 5 = f(x) es una parábola convexa. Sus puntos de 
corte con el eje OX son: 


=l 


2 = 
21 — Tr +5=0 e i gs 


y = —x? + 8x — 7 = g(x) es una parábola cóncava. Sus puntos de 


2 
=2Xx^-7X+5 
corte con el eje OX son: i 


4 


Puntos de corte de las dos parábolas: 


2x? —Tr+5 = —r?+8r—7 e 3z? —15r+12=0 e i a 


En consecuencia, al área que se pide será 
4 
A= | (90) $0) de 
1 


Calculamos una primitiva de g(x) — f(x): 


fut -Joj de=- J (32? — 152 +12) de =- (+ 2 122) | 


luego: 


Luego, finalmente, 
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31. Calcular el área encerrada por la gráfica de la función y = z? +1, la recta horizontal y = 1 y la recta vertical 
DO 


La curva y = z? + 1 corta a la recta y = 1 en un sólo punto: 
(x,y) = (0, 1). 


Además, en el intervalo [0, 3], la curva y = z? +1 está por encima 
de la recta y = 1 (x? +1 > 1 VY € [0, 3)). 


En consecuencia, el área que se desea calcular es: 


3 3 
isl (0 +11) a= f a? dz. 
0 0 


Una primitiva de x* es 
1 
3 4 
x? dr = -x 
J 4 


Por lo tanto, 


1 
32. Calcula el área de la región limitada por la curva y = x e7? y las rectas x = 3 ey=0 


La función f(x) = xe%% es positiva en (0,—00), negativa en 


(oo, 0) y se anula en z = 0. y 
Luego la región delimitada es la mostrada en la figura y su área 


viene dada por: 
1/3 
A= J xe 5% da 
0 


Calculamos en primer lugar una primitiva del integrando: 


F(z) = I ze da 


Esta integral indefinida se calcula por partes: 


v =e ; U=- 7e 
1 1 1 1 
=e Wma Tes n~” (e+ 5) 


Una vez calculada la primitiva, utilizamos la Regla de Barrow para calcular el valor de la integral definida: 


30 
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33. Calcular el área de la región del plano encerrada entre las gráficas de las curvas: 
g= mw ntl, yor +rord 


y las rectas verticales x = 1 y x = 5 


Denotaremos f(x) =2xz?—-x+1 y glæ)=r?+r+4 

Se trata de dos parábolas convexas, ya que el coeficiente principal es, en 
ambos casos, positivo. 

Para ver dónde se cortan planteamos la ecuación: 


—1 
22° —-z+1=2+r+4 > 122-21-3=0 e a= 3 
luego, en el intervalo [1,5], las dos parábolas se cortan en el punto de 
abscisa x = 3. 
Es fácil ver que la situación es la que se refleja en la figura adjunta. En 
consecuencia, 


-1 J 1 3 5 


3 5 
A= f (gla) — H(x)) de + f (Ha) — g(2)) de = A1 + A? 


Calculamos una primitiva de f(x) — g(x): 


F)= [o -gade = f (2 TEP d)d2= f (è - 22-3) dr =È -2° -32 


(Obviamente, una primitiva de g(x) — f(x) es, entonces, —F(x)). 


Ai -f (glz) — f(x) de = -ro a a 
Az = / (Hu) — glz)) dz = Fo]; T 
a244 6,328 

1 09 = 33. 3 


34. Calcular el área de la región del plano encerrada entre las gráficas de las curvas: 
y=a =m y=l=y 


y las rectas verticales x = 0 y x = 2 


— == 
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Denotaremos f(x) =1?—x, y g(a)=1-—z. 


Se trata, respectivamente, de una parábola convexa (ya que su coeficiente 
principal es positivo) y de una recta. Para ver dónde se cortan planteamos 
la ecuación f(x) = g(x), es decir, 


-1 


ae=1-% e 142-1=0 a 1 


Luego, en el intervalo [0, 2] ambas curvas sólo se cortan en el punto de 
abscisa x = 1, y es fácil ver que las gráficas son como en la figura adjunta. 
En consecuencia se tiene: 


A= | (0) Hadr+ | 9D) de = A +4 


Comenzamos por obtener una primitiva de f(x) — g(x): 


F)= fuga | (è -I)a =$ -2 


Obviamente, una primitiva de g(x) — f(x) es, entonces, —F (x). Se tiene, pues: 


“= | (ato 


x 
| 
= 
ea, 
3 
x 
a 
3 
Il 
ss | 
| 
3 
3 
is 
Il 
| 


35. Calcula razonadamente el área de la figura plana limitada por las curvas 
Ha == 20 ES y gla) =a” = 2 = 13 


y las rectas verticales z = —3 y x = 1. 


La gráfica de f(x) es una parábola cóncava (ya que su coeficiente prin- 
cipal es negativo) y la de g(x), por el contrario es una parábola convexa. 
Para ver dónde se cortan planteamos la ecuación f(x) = g(x): 


124 204+3=124?-20-13 & 22-22-8=0 > T. E 
Luego, en el intervalo [—3, 1] ambas curvas sólo se cortan en el punto de 
abscisa x = —2, y es fácil ver que las gráficas son las de la figura. En 
consecuencia se tiene: 


a=. (alo) Fa) de+ | ($) gla) de = A + Az 


Comenzamos por obtener una primitiva de f(x) — g(x): 
—2 
F(x) = fiw — g(x))dx = J (—2x? + 4x + 16) de = Ta + 22? + 167 


Obviamente, una primitiva de g(x) — f(x) es, entonces, —F (x). Se tiene, pues: 


e 


Matemáticas Aplicadas a la Biología Colección de ejercicios resueltos. Tema 3 


Ke a -f@)da = |-Fe)] =F y 4&= I Ñ fa) -glade = [Fo] =30 
A= +A = +36= a 


36. Dada la función f(x) = 1?(e* — 1) 
a) Estudiar su signo. 
b) Calcular una primitiva. 


c) Calcular, razonadamente, el área de la región del plano limitada por la curva y = f(x), el eje OX y las 
rectas verticales x = —1 y x = 1. 


a) Comenzamos por estudiar los puntos de corte de la gráfica de la función con el eje OX. Para ello planteamos 
la ecuación: 
r =0 & r=0 
r’ (e —1)=0 & z=4 o bien 
er-1=0%3 2=0 


Luego el único punto de corte es x = 0. Puesto que x? > 0 Vx e R, 


fa)<0 > e -1<0> e<1l > x<0 z f(x) < 0 en (—oo, 0) 
f(£)>0 & e-1>0 >l zr>0 F 


b) Calculamos una primitiva de f(x). Integramos por partes: 


v=e"—-1; v= -r 


F@)= | 1oar= f èE- 1d = u= w=2 [ea fajas 


La nueva integral que aparece también se calcula por partes: 


u = T; u'=1 q? y? 
frutero) de =2 | e(a) dx = at | =2|x| e” - — -f e — — 
v=e=x, v=- 2 2 


Luego, sustituyendo arriba se tiene: 


x? z3 gs 
F(x) = [trar era) 2x G 3 ) 2e” y 3 | = g'e”? —2re" +a? +2e"— = 


z Pa 2042) - 


Matemáticas Aplicadas a la Biología Colección de ejercicios resueltos. Tema 3 


c) Del apartado a) se deduce que la región cuya área hay que calcular 
es la de la figura adjunta, y que por lo tanto su área es: 


1 


as I fa) a f Ha) de = — [ræ] +r) ¿= 4+4 
0 


A, =-|F(2)]_. =—[F(0) - F(-1)] = F(-1) - F(0) = 


-1 


e 


=M 


Ejercicios del Tema 4: Métodos numéricos 


37. La medición de la temperatura de una habitación en un día proporcionó la siguiente tabla de valores: 


Tiempo (horas) [E 2 8 
Temperatura (Celsius) | y | 12 18 24 23 


Aproxima mediante interpolación lineal a trozos los datos anteriores, calculando la aproximación de la 
temperatura de la habitación a las 4 horas 15 minutos. 


Puesto que la variable independiente x se mide en horas, tenemos que pasar el dato 4 horas 15 minutos a horas: 
x = 4.25 horas (15 minutos es la cuarta parte de una hora). 
La interpolación lineal a trozos consiste en interpolar mediante la función polinómica obtenida al unir cada dos 
puntos consecutivos mediante una linea recta. 
Para calcular el valor de esta función en x = 4.25 hay que comenzar determinando entre qué dos valores de x 
se encuentra: en este caso entre x = 3 y x = 5. 

38 


25 
28 
15 


18 


=i 8 1 2 3 4 5 6 


Ahora hay que calcular la ecuación de la recta que pasa por los puntos (3,24) y (5,23): 


23 — 24 —1 -3 -3 
y E A E ii 
5-3 2 2 
La aproximación buscada es el valor de esta función en z = 4.25: 
4.25-—3 1.25 
y (4.25) = 24 3 = 24 3 23.375 


38. Tras ser administrado, se ha medido la concentración en plasma de un fármaco en intervalos de cuatro horas, 
obteniéndose: 
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4 8 12 16 20 24 


Tiempo ¿|0 
A a 


t 
Concentración | C 


Calcula una aproximación (mediante el método de los trapecios) de la cantidad total de fármaco absorbido 


por el organismo, es decir, 
24 


C(t)dt. 
0 


La fórmula de los trapecios para aproximar una integral definida es: 
b b—a 7 
FO de E [F(4) +29) Alta) + Flaa) 
E de k=2 


donde t1, ta, ...t, 1 son n +1 puntos regularmente espaciados del intervalo [a, b]. 
En este caso el intervalo es [0,24] y n = 6, ya que los 7 puntos dados forman 6 subintervalos de igual longitud 
h=4: 
Tiempo t ti ta t3 ta t5 te t7 
0 4 8 12 16 20 24 
Concentración | C | 0 18 32 23 10 4 0 


Luego se tiene: 


2 6 
J Op dir e [O(t1) +29 Clt) + Cltnsa)] 
0 


x6 
k=2 


2 ” [co + 2(0(4) + C(8) + C(12) + C(16) + 0(20)) 7 c(2a) 


=2|0+2(18+32+23+10+4) +0 = 4 x 87 = | 348 


39. La medición de la temperatura de una habitación en un día proporcionó la siguiente tabla de valores: 


Tiempo (horas) |z| 0 2 4 6 
Temperatura (°C) | y | 16 19 25 21 


a) Utilizando interpolación lineal a trozos, calcula una aproximación de la temperatura de la habitación 
a las 3:30 horas. 

b) Calcula una aproximación de la temperatura media (promedio integral) de la habitación en el intervalo 
de tiempo [0, 6], utilizando la fórmula de los trapecios. 


a) Puesto que la variable independiente x se mide en horas, tenemos que pasar el dato 3 horas 30 minutos a 
horas: x = 3.50 horas. 
La interpolación lineal a trozos consiste en interpolar mediante la función polinómica obtenida al unir 
cada dos puntos consecutivos mediante una linea recta. 
Para calcular el valor de esta función en x = 3.50 hay que comenzar determinando entre qué dos valores 
de x se encuentra: en este caso entre x = 2 y x = 4. 
Ahora hay que determinar la ecuación de la recta que pasa por los puntos (2,19) y (4,25): 


25-19 


=19 
Y * 43 


(x — 2) = 19 + 3(x — 2) 


= 20 — 
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La aproximación que buscamos es el valor de esta función en x = 3.50: 
y = 19 + 3(3.50 — 2) = 19 + 3 x 1.50 = 19 + 4.50 = | 23.50 
b) El promedio integral de una función f(x) en un intervalo [a, b] se define como: 


1 b 


En este caso se quiere calcular 


Aproximamos el valor de la integral en la fórmula anterior mediante la fórmula de los trapecios 
(la, b] = [0, 6], n = 3): 


6 
6—0 
0 f(æ)de x 5 (16 p2x 194 2x 25421) = 12% 


de donde 


- 1% 
F ~ = = [20.8333 


40. La evolución en tiempo del número de individuos de cierta especie (medido en miles) sigue la ley P, (t) = t, 
donde t es el tiempo (medido en años). La población de una especie competidora sigue la ley P2(t) =2—e?. 
Se considera que ambas especies tienen su supervivencia asegurada cuando coinciden en número. 


a) Mostrar razonadamente que el instante en que las especies coinciden en número tiene lugar entre t = 0 
y i=l, 


b) Aplicar el método de bisección para calcular aproximadamente el instante de coincidencia con un error 
menor que una décima. 


a) El instante en que ambas especies coinciden corresponde al valor de t para el que se tiene 
P(t) = P(t) => t=2-e = f(t) =t-2+e=0 
Se tiene: 
F(0) =-2+1=-1<0 y FU) =1-24+e=e-1>0 
lo que implica, por el Teorema de Bolzano, que f tiene al menos un cero en el intervalo [0, 1]. 


Además, puesto que f'(t) = 1+e? > 0 para todo t € R, la función f es monótona creciente, y en 
consecuencia, no puede tener más de un cero en dicho intervalo. 


En resumen: | existe un único punto en el intervalo [0, 1] en el que coinciden P, y Pa|. 


b) Aplicamos el método de bisección en el intervalo [0,1] para aproximar el cero de la función f(t): 


a b fa) FO) tm £ (tm) 


0 1 -1 1.7183 0.5 0.1487 

0 0.5 0.1487 0.25 -0.4660 
0.25 0.5 -0.4660 0.375 -0.1700 
0.375 0.5 -0.1700 0.4375 


= 37 — 
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El último intervalo, [0.375, 0.5], tiene longitud= 0.5 — 0.375 = 0.125. Luego, dando como aproximación su 
5 
= 0.0625 < 0.1. 


punto medio, t3 = 0.4375, cometemos un error menor que 


Podríamos haber decidido a priori el número n de iteraciones a realizar utilizando la fórmula que nos dice 
que para conseguir un error menor que e, debe ser: 


b— 
(24) 
—— = -1 
>) 
que, en este caso indica: 
E mm (3) — In(10) De 
"- Ta) mao 225 


que confirma lo dicho antes. 


41. Se desea calcular 


2 
| x lng dgr 
1 


a) Calcula una aproximación del valor de esta integral mediante la fórmula de los rectángulos (extremo 
izquierdo) utilizando 4 subintervalos. 


b) Calcula también una aproximación mediante la fórmula de los trapecios, asimismo con 4 subintervalos. 


1 
Dividimos el intervalo [1,2] en 4 subintervalos de igual longitud h = I~ 0.25: 


Ti T2 T3 T4 T5 
1 2 


Calculamos los puntos y los valores de f en ellos: 


k č Tk (xy) 
1 1 0 

2 1.25 0.2789 
3 1.5 0.6082 
4 1.75 0.9793 
5 2 1.3863 


a) La fórmula de los rectángulos: 


froar= f imadem (se) + tea) + fea) + ten) 


= 0.25(0 + 0.2789 + 0.6082 + 0.9793) 


= 0.4666 


b) La fórmula de los trapecios: 


frtyar= f zile) de = (101) +21(02) +24 les) +2/(04) + F(25)) 


1 
=z (0 +2 x 0.2789 + 2 x 0.6082 + 2 x 0.9793 + 1.3863) 


= | 0.6399 
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42. La medición de la temperatura de una habitación en un día proporcionó la siguiente tabla de valores: 


Tiempo (horas) |=| 0 2 4 6 
Temperatura (°C) | y | 16 19 25 21 


a) Utilizando interpolación lineal a trozos, calcula una aproximación de la temperatura de la habitación 
a las 3:30 horas. 


b) Calcula una aproximación de la temperatura media (promedio integral) de la habitación en el intervalo 
de tiempo [0, 6], utilizando la fórmula de los trapecios. 


a) La interpolante lineal a trozos de los datos contenidos en la tabla tendrá la gráfica reflejada en la figura. 


Para calcular su valor en un punto x hay que: 
(1) Identificar el intervalo [x;, £i+1] en el que está el punto. 
(2) Evaluar en x la función lineal que interpola los valores (2;,Y;), (UT;i+1, Yi+1)- 


En este caso, queremos evaluar la temperatura a las 3:30 horas, es decir, en x = 3.5. Este punto pertenece 
al intervalo [22,13] = [2, 4]. 
Calculamos ahora la ecuación de la recta que pasa por los puntos (2,19) y (4, 25): 

y-19 2 


Evaluamos y = 3x + 12 en z = 3.5: y = 3 x 3.5 + 12 = 22.5. 


Luego, | el valor estimado para la temperatura a las 3:30 horas es y = 22.5°C |. 


b) Por definición, el promedio integral de una función f(x) en un intervalo fa, b] se define como: 


_ 1 b 
5 | toa 


En este caso, [a,b] = [0,6] y b — a = 6, luego: 


=i TEGE 


Utilizando la fórmula de los trapecios con 3 subintervalos, podemos aproximar la integral entre 0 y 6 de 
la temperatura: 


6 


10) de m 3 (Flen) + 2f (22) + 2f (2) + F(24)) 


0 
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b-a 6 
=D, 
3 


Tenemos, por tanto: 


donde h = 


6 
/ f(x)dx =16+2x 19+2x25+21 
0 


= 16 + 38 + 50 + 21 = 125 


Luego, finalmente, tenemos, para la temperatura media en el intervalo de tiempo [0,6] el siguiente valor: 


=l 
F = z 125 ~ 20.83 


43. La evolución en tiempo del número de individuos de cierta especie (medido en miles) sigue la ley P, (t) = t, 
donde t es el tiempo (medido en años). La población de una especie competidora sigue la ley P2(t) =2—e?. 
Se considera que ambas especies tienen su supervivencia asegurada cuando coinciden en número. 


a) Mostrar razonadamente que el instante en que las especies coinciden en número tiene lugar entre t = 0 
yu= il 


b) Aplicar el método de bisección para calcular aproximadamente el instante de coincidencia con un error 
menor que una décima. 


a) Queremos justificar que existe algún punto entre t = 0 y t = 1, es decir, en el intervalo [0,1], en el que 
ambas funciones coinciden: 
P(t) = P(t), ie t=2—e!. 


O, lo que es lo mismo, queremos justificar que la ecuación 
t=2-e => t-2+e=0 


tiene alguna solución en el intervalo [0, 1]. 


Para ello utilizamos el Teorema de Bolzano. Denotamos f(t) =t—2+e!. Esta función verifica: 
F(0) =-2+1=-1<0, FU) =1-24+e=e-1>0. 
Es decir, f cambia de signo en [0, 1], en consecuencia, tiene al menos un cero en dicho intervalo. 


b) Utilizamos ahora el método de bisección para aproximar dicho cero: en cada iteración, dividimos en dos 
partes iguales el intervalo correspondiente, quedándonos con la parte en que se produce el cambio de signo. 
Para asegurar un error menor que una décima (0.1), detendremos este proceso cuendo tengamos el cero 
localizado en un intervalo de amplitud menor que 0.2, y tomaremos como aproximación el punto medio 
de ese intervalo. 


a |b| sa) | SO [| b-a| tu | Slem) | 
0 1 | —1<0 | 1.7183 >0 1 0.5 1.1487 > 0 
0 0.5 <0 >0 0.5 0.25 —0.466 < 0 
0.25 | 0.5 <0 >0 0.25 0.375 —0.17 < 0 
0.375 | 0.5 <0 >0 0.125 0.4375 
Observamos que el último intervalo [0.375, 0.5] tiene amplitud = 0.125 < 0.2. En consecuencia, su punto medio, 
A12 
tin = 0.4375 proporciona una aproximación del cero de f(t) con un error < a = 0.0625 < 0.1, como nos 


piden. 


—4A0 


Ejercicios del Tema 5: Ecuaciones 
diferenciales 


44. Para la ecuación diferencial 


se pide: 


a) Obtener la solución general. 


b) Determinar la solución que verifica y(3) = 2. 


a) Esta ecuación es de variables separables: 


2 
3 y +1 MY, 1 0 2y f 1 
= > = < —=— dy= | —— d 
x(x — 1) ytl” x(x — 1) y2? +1 sd a(x—1) PA 
La integral del segundo miembro se calcula escribiendo el integrando como suma de fracciones simples: 


1 A B Es 
+= 


2yy 


xa-1) z zs-—l B= 1 


En consecuencia 


I 1 n 
m +i= f (++) dz = —Inļe|+n|z — 1| +C = m|? 


1 
+C 
De aquí, tomando exponenciales en ambos miembros, se tiene 


x— 1 


ln +C ln 


y +1l=e 


— 1 —1 
ep =c(2 )-1 e 9-2 JE )- 


b) Veamos ahora para qué valor de C se tiene y(3) = 2. Lo primero que se observa es que, puesto que 2 es 
positivo, hay que tomar el signo positivo de la raíz cuadrada: 


T 


— — 1 
2= 40? -1-42 2 > E 3 s 1243=2C => o=" 


Así pues, la solución buscada es: 


41 
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45. Para la ecuación diferencial 
2yy' = (y? + 1)x cos(x), 


se pide: 


a) Obtener la solución general. 


b) Determinar la solución que verifica y(0) = 1. 


a) Esta ecuación es de variables separables: 


2y 2y 
c u2 = = 
2yy' = (yf + 1)z cos(x) > +1 y =xcos(x) > i P4i dy = J x cos(x) dx 


La integral del segundo miembro se calcula utilizando la fórmula de integración por partes: 


frost) dx = | la H e = g sen(x) — J sen(x) dz = x sen(x) + cos(x) + C 


En consecuencia se tiene 
In |y? +1| = x sen(x) + cos(x) + C 
De aquí, tomando exponenciales en ambos miembros, se tiene 


y? +1= e sen(x)+cos(1) +0 _ e? sen(x)+cos(x) eC = Ce? sen(x)+cos(x) 


s y? = Ce sen(x)+cos(x) _ lo y= YC er sen(x)+cos(2) — 1 


b) Veamos ahora para qué valor de C se tiene y(0) = 1. Lo primero que se observa es que, puesto que 1 es 
positivo, hay que tomar el signo positivo de la raíz cuadrada: 


2 


1 = +y C e0sení0)+cos(0) — 1 = +yCe—1 => 1=Ce-1 > 2=Ce e C 


Así pues, la solución buscada es: 


y=+ i ex sen(s)+cos(x) — 1 
e 


46. Resuelve el siguiente problema de valores iniciales 


y 

a E 

y o 
y(0) =1 


i y 1 1 
= _ dy = de > Inly|=ln ja +1]|- = 1 
y cd e f dy J= z n|y| =1n lx I4C => yn =Clx+1) 


Para obtener una solución particular de la ecuación dada, calculamos ahora 


K= [oge fee f (1-5) dezhail 
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de modo que la solución particular buscada es 
ylz) = K(0)G(2) = (2 -mle +1)) (2+1) 
y la solución general de la ecuación lineal es, finalmente, 


y=C(2+1) + (2 Inje+1)) (1 +1) = (C+z-m|e+1]) (£x +1) 


Para hallar la solución del problema de valor inicial, imponemos la condición inicial: 


1=y(0)=(C t0- In]0 4 1) (0+1)=C-In(1)=C 
Luego la solución pedida es: 


y=(1+2-Inje+1)) (1 +1) 


47. Resuelve el siguiente problema de valores iniciales 


Se trata de una ecuación lineal. Calculamos en primer lugar la solución general de la ec. homogénea asociada: 


R Y 1 1 
= > |-dy= d l =l 2| +C < =C 2 
y= 4 0 [iws | i e nisme m=C(+2) 


Para obtener una solución particular de la ecuación dada, calculamos ahora 


1 1 2 1 
K(x) = | b(2)— dz = —— dr = 1— =x-—2l 2| = ln 5 
(x) I lag” pe H 3) dx =x n|e+2| =+ E 


de modo que la solución particular buscada es 


yp(x) = K(2)G(x) = (z —2ln|e + 21) (z +2) = (e +In ==) (£ +2) 


y la solución general de la ecuación lineal es, finalmente, 


y = Clx +2) - € 21n [2 + 2)) (1 +2) = (0 +2 -21n|2+2)) (x +2) = (c Hg +ln : 2) (z +2) 


Para hallar la solución del problema de valor inicial, imponemos la condición inicial: 
l1=y-1)=C-1 & C=2 


Luego la solución pedida es: 


y=(2+2 mje +2) @+2)= (240410 z )@+2 
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Matemáticas Aplicadas a la Biología Colección de ejercicios resueltos. Tema 5 


48. Resolver el siguiente problema de valores iniciales: 


f y' +2 tg(t)y = sent 
y(0)=4 


1 
y =-2 tewe | dy =-2 f tgtdt o im y = 21n|cost| + C = ln cos? t + C > yn(t) = C cos? t 
y 


1 sen t 1 1 
a K= J sent dt = dt = — t —— cos? t = cost 
J cos? t J cos? t cost Yplt) cost 


= | y(t) = C cos? t + cost 


= 4 = y(0) = C cos? 0 +cos0 =C +1 => C=3. Luego la solución es |y = 3 cos?t+ cost 


49. a) Calcular la solución del siguiente problema de valor inicial: 


y +y= 
y(0) =0. 


b) Estudiar el comportamiento de la solución del problema anterior cuando t > +00. 


Se trata de una ecuación lineal: . 7 


e 
é y =-y+ 


li = 
Yy TY 4-2 


Ec. homogénea: 
1 
y = —y > Pu=- fa <> lay = -t+ C 5S Yn = C e™ 
Y 


Solución particular de la completa: yp = K(t)e”* con 


eto 1 1 
(+) las i Pa: 


Para calcular la integral hay que escribir la fracción como suma de fracciones simples. Para ello, se comienza por 
factorizar el denominador: 4 — t? = (2 — t)(2 + t). Puesto que se descompone en dos factores simples distintos, 
la descomposición que necesitamos será de la forma: 


1 1 A B 


¿2 (2-500Q+t) A 


Multiplicamos en ambos miembros por 4 — t? o, lo que es lo mismo, por (2 — t)(2 + t), y obtenemos: 
1=A(Q+t)+B(2-t) 


Necesitamos encontrar A y B para que la anterior igualdad sea cierta Yt € R. En particular dando valores a t 
se tiene: 


t=2 1=4A A=1/4 
t=-2 1=4B B=1/4 


Luego, se tiene finalmente, para la integral: 


1 1/4 1/4 1 1, 12+t 
dis ( ) at = ( n|2—+| +10 12 8) =31 a 
[r l 2-1 Da a MAA m= 


=j 
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Luego la solución particular buscada es:: 


2+t 


=K(t = in a 
Yp (He~ 2: e 
y la solución general de la ecuación diferencial es: 
1, ¡2+t 
y=" Het oen 


Esto es la solución general de la ecuación diferencial. Debemos ahora calcular la solución particular que verifica 


la condición inicial dada: 
1 2+0] o 
y(0)=0 & 0= 1” 30 


Luego la solución del problema de valor inicial dado es: 


l Ce = W(1):1+0-1=0 


Por último, el comportamiento de la solución cuando t > +00 se estudia calculando su límite cuando t => +00. 
Se tiene 


2+t 
E =1-11=1> lm In | =0 
Means 2 t>+o0 

y 

lím e*=0 

t>+00 
En consecuencia 

1 2+t 

lím -ln Teje =0 
t>+00 4 2-t 


50. Debido a un vertido tóxico, las aguas de un lago, que inicialmente estaban limpias, comienzan a ser conta- 
minadas. Sea y(t) la cantidad de contaminante (gramos) por m? en el instante de tiempo t (medido en días). 
Se sabe que y(t) sigue la ley 


1 
1 
= — k = 
y =¿(k=y) 
donde k es una constante desconocida. Pasados 4 días tras el vertido, la contaminación es de 86.5 gr/m. 


a) Calcular la función y(t) y el valor de la constante k. 


b) Calcular la cantidad de contaminante por m? que tendrán las aguas del lago 10 días después del vertido. 


a) La solución general de la ecuación diferencial es y = k — Ce™*/?: 


1 1 
——dy= | -dt > —In|k — y| = 
Fw Jus n|k — yl 


Puesto que, inicialmente, no había contaminación, debe ser: 0=y(0)=k-C => C= 


1| t 
"ky 


Puesto que 4 días más tarde era y = 86.5, debe ser: 


LLR 
86.5 = y(4) = k — ke7? = k(1- e?) => k=C= T 


= 100. 
Así pues la función que da la concentración de contaminante t días después del vertido es: 


y(t) = 100(1 — e7/2) 


b) Conocida y(t), se puede calcular la cantidad de contaminante que habrá en el lago 10 días después del 
vertido: 
y(10) = 100(1 — e7 1/2) = 100(1 — e75) ~ 99.32 


— 45 
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51. Una determinada especie ha sufrido una epidemia, con lo que la tasa de mortalidad ha superado la tasa de 
natalidad. Se sabe que y(t), el número de individuos vivos de la especie en el instante de tiempo t, es la 
solución del siguiente problema de Valor Inicial: 


donde a > 0 es un parámetro e yo > 0 es el número de individuos de la especia presentes al comienzo del 
periodo de observación. 


a) Determinar la solución al problema anterior para a = 1, yo = 10. 


b) Probar que, para tiempos grandes, la especie tiende a la extinción. 


a) Debemos resolver la ecuación diferencial: y’ = —y—y? = —y(1—y) que claramente tiene las soluciones 
triviales constantes y = 0 e y = —1 


Se trata de una ecuación logística, que re-escribimos de la forma (para y # 0 e y # —1): 


dy dy dy 

~= =-y(l +y) e = -dt i z fa 

dt y( +y) y(1 +y) 
1 1 

Descomponemos z= como suma de fracciones simples, de la forma : 
+4? yl+y) 
1 A B y= > 1=A 
= S1=A(1+y)+B 
y(l+y) y 1+y Ey di 1 -1 > 1=-B 


Se tiene, pues 


E aF 
= dy=— | d =  Inly|—In|1+yl|= ln 
Jna y 1+y W= iatt 1 


Despejamos la y: 


Ta =Ce* e y=Ce*(1+y=Ce*+Ce "y e (1-Ce™)y=0e™ 
y 
luego 
= o, lo que es lo mismo Y= pa] CER 
Sustituyendo y(0) = 10, obtenemos: 
1 1 1 1 1+10 11 
= 10 = y(0) = - > C-1l=— C H1= = 
dé WOS apin g-i 10 10 10 10 
luego la solución del problema de valor inicial es: 
DEI E: + 
"= dé=1 UD, , me-10 
—et-1 
10 


b) El comportamiento de y(t) para tiempos grandes se obtiene calculando el siguiente límite: 


í ()= í 10 E 
o l1et -10 


0 


lo que significa que la especie tiende a la extinción. 
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52. Las aguas de un estanque de regadío se encuentran contaminadas, debido a los plaguicidas. El estanque 
contiene 3000 m? de agua con una concentración de 2 g/m? de DDT. Para purificar el estanque, se comienza 
a introducir agua limpia a razón de 42 m*/h. Simultáneamente, el agua, que se supone mezclada de forma 
instantánea, se vacía a la misma velocidad. 


Para controlar la operación, se necesita conocer la función que proporciona la cantidad de DDT que hay 
presente en el estanque en cada instante del proceso. 


a) Determina la ecuación diferencial que verifica dicha función. Determina también la condición inicial 


que hay que imponer para calcularla. 


b) Resuelve el problema de valor inicial planteado en el apartado anterior. 


c) ¿Cuánto tiempo debe pasar para que la concentración de DDT en el estanque descienda hasta 0.5 g/m3? 


Utilizaremos las siguientes unidades: 
Volumen: metros cúbicos (m3) 
Cantidad de contaminante: gramos (g) 
Tiempo: horas (h) 


a) 


b) 


Puesto que entra agua al estanque a la misma velocidad a la que sale, en el estanque hay siempre el mismo 
volumen de agua: 
V = 3000 m? 


Si y(t) es la cantidad total de DDT que hay en el estanque en el instante t, entonces, puesto que y’ es la 
variación de y por unidad de tiempo, se tendrá: 


y'(t) = cantidad que entra/hora — cantidad que sale/hora 


Cantidad que entra por hora: el agua que entra en el estanque está limpia. Por tanto, entran en el 
estanque 0 g de DDT por hora. 


Concentración de contaminante en el estanque en el instante t: será igual a la cantidad total 
de DDT que hay en el estanque en el instante t, es decir, y(t), dividido por el volumen total de agua que 
hay en el estanque en el instante t, que en este caso se mantiene constante e igual a 3000 m*, como hemos 
y(t) 
3000 
Cantidad que sale por día: salen igualmente 42 m*/h, con una concentración de contaminante igual, 


42 
por lo dicho antes, a Y g/ m. En consecuencia, salen de la piscina - 


Por lo tanto, ya se puede establecer la ecuación diferencial que verifica y(t): 


visto antes. Es decir, la concentración en el instante t es g por cada më. 


g por hora. 


4y 4 


E 
y =0= 3000 T000 


La condición inical será la que indique la cantidad inicial de DDT en el estanque. Puesto que inicialmente 
hay 3000 m*de agua con una concentración de 2 g/m, la cantidad inicial será 


y(0) = yo = 2 x 3000 = 6000 (g) 
Resolvamos, en primer lugar, la anterior ecuación diferencial, que es de variables separables: 


14 1 14 14 
'=-— ye | “d=-"— | dt > Inly =- —t+C & y= C e714/1000 
E 1000 ” e id el nlyl == Toot Y ARE 


La condición inicial: 
y(0) = 6000 


EE p an 
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Sustituyendo en la expresión de la solución general obtendremos el valor adecuado de C y, con éste, la 
solución particular que verifica la condición inicial: 


6000 = y(0) = C 


Luego la función buscada es 


y = 6000 e” 14/1000 


c) La función y(t) calculada en el apartado anterior nos da la cantidad total de DDT que queda en el 
estanque en el instante t. Puesto que el volumen total de agua en el estanque es 3000 më, tenemos que la 
concentración de DDT en el estanque en el instante t es: 


Tenemos que calcular el valor de t para el que se tiene 


c(t) =0.5 
es decir, 
0.5 14t 
0.5 =2 —14t/1000 == —14t/1000 a 0 25 E e 1 0 5 
E e 2 1000 ~ (025) 
1000 


53. Se considera una población que evoluciona según una ecuación logística pero que, además, está sometida a una 
depredación que reduce el tamaño de la población en 10 individuos en cada tiempo t. Este comportamiento 
hace que la velocidad de cambio de la población se exprese de la siguiente forma: 


O E, 
y = yl 40) ; 


donde y(t) representa el número de individuos en el tiempo t, y(1 — y/40) expresa el término logístico y el 
segundo término expresa la depredación. 


a) Si inicialmente hay 100 individuos, ¿cuántos individuos habrá en cada tiempo t? 


b) ¿Se extinguirá la población? 


a) El problema de valor inicial que hay que resolver es: 


Calculamos en primer lugar la solución general de la ecuación. Comenzamos por factorizar el segundo 


miembro: 
Y 1 1/2 Z-t 2 
y(1— E) - 10 = zg (y(40— y) — 400) = 77 (1? — 40y + 400) = z5 (v - 20) 
Luego 
Mia 2? = |! J dt 
=p (y — 202” 40 
—1 —t 1 t t+C 
<= = C <= = PO = 
y-20 m" y=20 40 40 
40 40 


—A8g— 
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Imponemos ahora la condición inicial 
40 40 40 1 
y(0) 00 0+ C 80 o C 30 C 3 
Luego la solución del problema de valor incial es: 
40 80 

=20 = =P 

il pe +3 THp 

p nae T . 80 
El número de individuos de esta población en el instante t es y = 20 + AF] 
80 
2t+ 1 


b) Para analizar si la población se extingue observamos en primer lugar que la función y = 20 + 2 
es decreciente, lo que significa que la población disminuye con el paso del tiempo. Calculemos su límite 


cuando t — o0: 80 80 


Es decir, la población disminuye de tamaño pero no se extingue, ya que el número de individuos tiende 


asintóticamente a 20. 


54. El Bario 133 (199Ba) es una sustancia radioactiva y se desintegra a una velocidad proporcional a la cantidad 
presente. La vida media (tiempo que tarda en desintegrarse la mitad de cualquier cantidad inicial) es de 


10.53 años. 
a) Determina la ecuación diferencial que modeliza la evolución temporal de la cantidad de Ba. 


b) Calcula la constante de desintegración del Ba. 
c) Averigua el tiempo que tarda 1 kg de 133Ba en reducirse a 1 g. 


a) La ecuación diferencial que modela la desintegración radiactiva y expresa que la desintegración se produce 


a una velocidad proporcional a la cantidad existente es: 
y =-Ay 


donde y(t) es la cantidad de !33Ba presente en el instante t (tiempo, medido en años) y A > 0 es la 
constante de desintegración o decaimiento del +9%Ba. Las soluciones de esta ecuación son de la forma: 


y=Ce* CER 


b) La constante de desintegración se calcula a partir de la vida media. Supongamos que inicialmente hay una 
cantidad yo de Ba: y(0) = yo. Puesto que la vida media es Vm es el tiempo que tarda en desintegrarse 
la mitad, se tendrá y(Vm) = 2% De estas dos ecuaciones se puede despejar A en función de la vida media: 

y(0) =y=Ce”** <= C=yo, estoes y=ype ?** 
1 In(2 
m EE A AN 
Vm 


Yo — 
z= = Ye AVm ¿> `= 
2 
Esta es la relación entre la vida media y la constante de desintegración de cualquier sustancia radiactiva 


En este caso, pues, se tiene que la constante de desintegración del 133Ba es 


(2) (0.0658 
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c) De los apartados anteriores se tiene que la función que nos da la cantidad de +9Ba presente en el instante 
t, si se ha comenzado con una cantidad yy es: 


aa —0.0658 t 
Yy=Y0€ 
Queremos calcular para qué valor de t se tiene 


1 1 
y(t) = 10000700658 =1 ¿=> e700658t = —— <= —0.0658t = ln (—) 


1000 
z 1 
asis ni 

t= 068 2 (000) 


=1|105 años 


